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Abstract 



Extending our previous results for artificial viscosity systems, we 
show, under suitable spectral hypotheses, that shock wave solutions of 
compressible Navier-Stokes (cNS) and magnetohydrodynamics (MHD) 
equations undergo Hopf bifurcation to nearby time-periodic solutions. 
The main new difficulty associated with physical viscosity and the cor- 
responding absence of parabolic smoothing is the need to show that 
the difference between nonlinear and linearized solution operators is 
quadratically small in H s for data in H s . We accomplish this by a 
novel energy estimate carried out in Lagrangian coordinates; interest- 
ingly, this estimate is false in Eulerian coordinates. At the same time, 
we greatly sharpen and simplify the analysis of the previous work. 
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1 Introduction 

A well-known phenomenon in combustion is the appearance of "galloping", 
"spinning" , and "cellular" instabilities of traveling detonation fronts, appar- 
ently corresponding with Hopf bifurcation of the background solution; see, 
e.g., |KSl LyZl , LyZ2 \TZ1\ |TZ2j and references therein. Mathematically, 



there are two distinct aspects of this phenomenon. The first is to verify the 
spectral scenario associated with Hopf bifurcation, consisting of a conjugate 
pair X±(e) = 7(e) + ir(e) of complex eigenvalues of the linearized operator 
about the wave crossing the imaginary axis from stable (negative real part) to 
unstable side as bifurcation parameter e varies from negative to positive val- 
ues. Here, e measures variation in physical parameters such as heat release, 
rate of reaction, or strength of the detonation. 

The second issue is to show that this spectral scenario indeed corresponds 
at the nonlinear level to Hopf bifurcation: i.e., apparition of nearby time- 
periodic solutions u a of approximate period 2tt/t(e), e = e(a), branching 
from the steady solution at e(0) = 0. What makes this nontrivial is the 
absence of a spectral gap between zero and the essential spectrum of the lin- 
earized operator about the wave, which circumstance prevents the application 
of standard bifurcation theorems for PDE, as in, e.g., [CJ IHe} IMM[ IVTj . 
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The first issue has been studied for detonations by numerical |KSj . for- 
mal asymptotic |Erl \FW\ IBMRj . and Evans function |LyZl , LyZ2| tech- 



niques. The second issue has been studied in a much more general setting in 
[TZ1[ ITZ2| ISSj . In particular, it was shown in |TZ2] for smooth shock solu- 
tions of conservation laws with artificial (strictly parabolic) viscosity that a 
spectral Hopf scenario implies actual nonlinear Hopf bifurcation. This result 
was recently sharpened by Sandstede and Scheel [SB] to include also expo- 
nential localization of the perturbed solution (the result of |TZ2] asserts only 
(1 + decay) and exchange of spectral stability. The results of [SSj were 

obtained by quite different "spatial dynamics" techniques, viewing the prob- 
lem as a spatial elliptic boundary value problem on the space of time-periodic 
functions, whereas the original results of [TZ2] were obtained by "temporal 
dynamics" techniques centered around the temporal evolution operator. 

Our purpose in the present paper is to extend the latter, nonlinear, re- 
sults to the physical (nonstrictly parabolic) systems of compressible gas- and 
magnetohydrodynamics (MHD) with either ideal or "real", van der Waals- 
type, equation of state by a modification of the original, temporal dynamics 
argument of |TZ2j . At the same time, we show that exponential localization 
may by an improvement of the central cancellation estimate be obtained by 
temporal as well as spatial techniques, recovering the sharpened estimates of 
[SSJ; moreover, this same improvement greatly simplifies the entire analysis. 
We expect, but have not verified that exchange of stability may be obtained 
by our techniques as well. 

In contrast to the detonation case, this investigation is not driven by 
known physical phenomena, but by the common mathematical structure with 
detonation. That is, in this case the mathematics suggests a possible phys- 
ical phenomenon, and not the reverse. It would be very interesting to look 
numerically or experimentally for its appearance, particularly in the rich set- 
tings of MHD or phase-transitional gas dynamics with van der Waals-type 
equation of state, where instability, similarly as in the detonation case (but 
in contrast to gas dynamics), frequently occurs [T]. 

Remark 1.1. A third approach to bifurcation in the absence of a spectral 
gap, besides those described in [TZ2\ ISSj . may be found in a recent work of 
Kunze and Schneider |KuS] in which they analyze pitchfork bifurcation in 
the absence of a spectral gap using weighted-norm methods like those used 
to study stability for dispersive systems |PWj . This approach has been used 
to treat stability of dispersive-diffusive scalar undercompressive shocks [Do]; 
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however, it does not appear to generalize to shock waves in the system case. 



1.1 Equations and assumptions 

Consider a one-parameter family of standing viscous shock solutions 

(1.1) U(x,t) = U £ (x), lim U £ (z) = U £ ± (constant for fixed e), 

z— >±oo 

of a smoothly-varying family of conservation laws 

(1.2) U t = T{e, U) := (B(e, U)U X ) X - F(e, U) x , U G M n , 
with associated linearized operators 

(1.3) L e := dF/dU\ u=f j* = d x B e (x)d x - d x A e (x) 

B £ (x) = B(U £ {x)), A £ (x) := Fu{U £ {x),s). Equations (Q are typically 
shifts B(e, U) = B(U), F(e, U) := f(U) - s(e)U of a single equation 

U t = (B(U)U X ) X - f(u) x 

written in coordinates x = x — s(e)t moving with traveling- wave solutions 
U(x,t) = U £ (x — s(e)t) of varying speeds s(e). Profiles U £ satisfy the 
standing-wave ODE 

(1.4) B(U)U' = F(e,U) -F(e,Ui). 
Denote 

(1.5) B £ ± := hm B £ (z) = F n (U±, e), A £ ± := lim A £ {z) = Fu(U e ± , e), 

Z—f±OC z— +±oo 

and A{e, U) = F u {e, U). 

We make the following structural assumptions. 

(Al) U = ( ^ ^ , B = (° Q ° b \ b nonsingular, where U G R n , 

U x G R n - r , U 2 G R r , and b G W xr ; moreover, the w-coordinate F^e, U) of F 
is linear in U (strong block structure). 

(A2) There exists a smooth, block-diagonal, positive definite matrix 
A°(e,U) such that A^An is symmetric and A® 2 b is positive definite but 
not necessarily symmetric (symmetric hyperbolic-parabolicity). 

To (Al)-(A2), we add the following more detailed hypotheses. Here and 
elsewhere, a(M) denotes the spectrum of a matrix or linear operator M. 
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(HO) F,B eC k , k>5. 

(HI) a(A £ u ) real, semisimple, nonzero, and constant multiplicity. 
(H2) cr(v4j_) real, simple, and nonzero. 

(H3) Rcr (i£ j dF(U±) - i 2 B{U±)) < —Mb, > 0, for f e R. 

(H4) Considered as connecting orbits of (11.41) . £7 e lie in an ^-dimensional 
manifold, £ > 1, of solutions (II. ip . obtained as a transversal intersection of 
the unstable manifold at U e _ and the stable manifold at Ul. 



Remark 1.2. Conditions (H1)-(H2) imply that U± are nonhyperbolic rest 
points of ODE (jl .4[) expressed in terms of the -^-coordinate, whence, by stan- 
dard ODE theory, 

(1.6) \di{U £ -U £ ± ){x)\ < Ce-^W, 0<£<k + l, 

for x ^ 0, some rj, C > 0; in particular, \(U £ )'(x)\ < Ce~ v ^. Condition 
(H4) implies in part that U e is either of standard Lax type, £ = 1, or non- 
classical overcompressive type, £ > 1, i.e., the hyperbolic convection matrices 
A± := Fjj{U±, e) at ±oo have, respectively, p — £ negative and n—p positive 
eigenvalues for 1 < £ < p < n. In the L clX CclSC, the (£ = l)-dimensional man- 
ifold of solutions (11 .ip consist of translates of U £ , and p is the characteristic 
family associated with the shock. For further discussion, see |ZH| IMaZT|IZ2j . 

Remark 1.3. Similarly as in [Z2, Z3J, it is readily checked that the conclu- 
sions of this paper hold also with (A1)-(A2) replaced by the weakened as- 
sumptions that there exist an invertible change of coordinates U — > (U\, W2) 
for which F\ is a linear function of U-y, W2 with du x Fi symmetrizable and 

B(U)d x U = (p^Q w ) with 9ttT((dW 2 /dU 2 )P) < -9 < and (HI) 

modified as in [MaZlJ to account for nonvanishing of B21, replacing An 
by A* := An — B2iB 22 1 Ai 2 ; see Remarks 13.51 and 1 5. 11 These assumptions are 
used to obtain the H s energy estimates of Section [3j 

Conditions (Al)-(A2) (or the alternatives described in Remark 13.51) are 
a slightly strengthened version of the corresponding hypotheses of |MaZ3t 
\Z2\ IZ3] for general systems with "real" , or partially parabolic viscosity, the 
difference lying in the strengthened block structure condition (Al): in par- 
ticular, the assumed linearity of the U\ equation. Conditions (H0)-(H4) are 
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the same as those in |MaZ3l |Z2| |Z3] . The class of equations satisfying our as- 
sumptions, though not complete, is sufficiently broad to include many models 
of physical interest, in particular compressible Navier-Stokes equations and 
the equations of compressible magnetohydrodynamics (MHD), expressed in 
Lagrangian coordinates, with either ideal or "real" van der Waals-type equa- 
tion of state; see Section \5\ for further discussion. The role of assumptions 
(Al)-(A2) in our analysis is discussed further in Appendix [A1 

A simple example is the equations of isentropic gas dynamics written in 
Lagrangian coordinates: 



;i-7) 



v t - u x = 0, 

Ut+Px = {(v/v)u x ) x , 



where v > denotes specific volume, u fluid velocity, p = p(v) pressure, and 
v > the coefficient of viscosity, with p v < corresponding to hyperbol- 
icity of the associated first-order system. Denoting Ui = v , U2 = u, and 
considering a family of traveling-wave solutions U(x,t) = U £ (x — s(e)t) with 
p v (v±) < 0, we obtain after the change of coordinates x — * x — s(e)t a family 
of stationary solutions U(x,t) = U £ (x) and equations 

/j g N f v t - s{e)v x - u x = 0, 

\ u t - s(e)u x +p x = {{v/v)u x ) x 

satisfying (A1)-(A2), (H0)-(H4). 

Remark 1.4. Equations (11. 7p written in Eulerian coordinates, or 



;i.9) 



Pt + {pu) x = 0, 

(pu) t + {pu 2 +p) x = uu x 



where p — 1/v denotes density does not satisfy (Al)-(A2), (H0)-(H4), vio- 
lating the second part of the strong block structure assumption (Al). (Here, 
U\ = p, U2 = m = pu is momentum, so that the first equation is still lin- 
ear; however, the parabolic term u xx = (m/p) xx is not in diagonal form.) 
Indeed, we make crucial use in the analysis of the Lagrangian formulation, 
as discussed in Appendix |A] 



1.2 Spectral criteria 

As discussed in |ZH[ IMaZl[ |ZT| |Z2| |Z3] , the linearized operators L £ have no 
spectral gap, since essential spectrum accumulates at A = on the imaginary 
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axis. In this situation, standard stability and bifurcation criteria based on 
isolated spectra of L £ are replaced by generalized versions expressed in terms 
of an associated Evans function. The Evans function D e (X), defined as a 
Wronskian of functions spanning the decaying manifolds of solutions of the 
eigenvalue equation 

(1.10) (L £ -\)u = 

associated with L £ at x —>■ +00 and x — ► — 00 is an analytic function with 
domain containing {9?A > 0}, whose zeroes away from the essential spectrum 
correspond in location and multiplicity with eigenvalues of L £ . Its behavior 
is also closely linked with that of the resolvent kernel of L £ , i.e., the Laplace 
transform with respect to time of the Green function G; see, e.g., |AG J\ IG"Z| 
IZH[ \ZS\ |ZT| |Z2] for history and further details. 

In [MaZlj there was established the following stability criterion. 

Proposition 1.5 ([MaZlj). Let U £ , (ll.2j) be a family of traveling-waves and 
systems satisfying assumptions (A1)-(A2), (HQ)-(H4). Then, U £ is linearly 
stable from L 1 — > L p if and only if there exist precisely £ zeroes (£ as in (H4)) 
of D £ (-) in the nonstable half-plane IRA > ; necessarily at the origin A = 0. 

We define an analogous Hopf bifurcation criterion as 

(T> e ) On a neighborhood of {3?A > 0} \ {0}, the only zeroes of D 
are (i) a zero of multiplicity £ at A = 0, and (ii) a crossing conjugate pair of 
zeroes A±(e) = 7(e) + zr(e) with 7(0) = 0, 0*7(0) > 0, and r(0) ^ 0. 

Remarks 1.6. 1. In the simplest, Lax case, under (Al)-(A2), (H0)-(H4), 
simplicity of the root A = 0, condition T> e (i), is equivalent to (ti £ , (u £ )') = 
ir e -(u £ + -u £ _) ^ for vr £ (constant) orthogonal to S(A £ _)UU(A £ + ) [GZllMaZ2l 
IZ1UZ2] . Under the normalization (ir £ , (u 6 )') = 1, operator U^f := (u £ )'(tt £ , f) 
plays the role of a "generalized spectral projection" onto 

KerL(e) = Span{(M £ )'}, 

and 7r e the role of a generalized left eigenfunction [MaZll IZH] . Note that n £ 
lies outside the domain of ITq, a consequence of the absence of a spectral gap. 
A similar, but more complicated condition holds in the overcompressive case 
[MaZTl lZTllZ2]. 

2. By Proposition Ol condition £> e with (A1)-(A2), (H0)-(H4) implies 
that U £ is linearly stable for e < and unstable for e > 0; that is, there is a 
transition from stability to instability at £ — 0. 
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1.3 Results 



We introduce the following notation, to be used throughout the paper. 

Definition 1.7. Let B 2 C B\ and X 2 C X\ denote the Banach spaces 
determined by norms \\U\\ Bl ■= \\U\\ H i, \\d x U\\ B2 ■= \\d x U\\ Bl + \\U\\ L i and 

\\U\\ Xl := \\e r ' {x) U\\ H 2, 
ll^C/IU.^II^t/H^ + lle^C/llifi, 

where r] > and (x) := (1 + Ix) 2 ) 1 / 2 . 

Our main result is the following theorem establishing Hopf bifurcation 
from the steady solution U e at e = under bifurcation criterion (T> e ). 

Theorem 1.8. Let U £ , (1 1.2ft be a family of traveling-waves and systems 
satisfying assumptions (A1)-(A2), (H0)-(H4), and (D e ). Then, for r > 0, 
f] > sufficiently small and C > sufficiently large, there is a C 1 function 
e(r), e(0) = ; and a C 1 family of time-periodic solutions XJ r (x,t) of (11.21) 
with e = e(r), of period T{r), T(-) G C 1 , T(0) = 27r/r(0), mft 

(1.12) C-V < ||U r - U e \\ Xl < Cr 

for all t. For Lax shocks, up to fixed translations in x, t, these are the only 
time-periodic solutions nearby in X\ with period T G [Tq,Ti] for any fixed 
< T < T\ < +oo; if 7^ U £ _, they are the only nearby solutions of the 
more general form XJ r (x — o~ r t, t) with XJ r (x, •) periodic. For overcompressive 
shocks, they are part of a C 1 {I — 1) -parameter family of solutions that are 
likewise unique up to translation in x, t. 

Remark 1.9. Bound (11.121) . by Sobolev embedding, includes also the result 
of exponential localization, |U r — U £ \ < Cre~ v ^. 

1.4 Analysis 

Theorem II .81 is proved using the general bifurcation framework established in 
|TZ2] together with pointwise Green function bounds established in |MaZl] 
for general hyperbolic-parabolic systems. This is for the most part straight- 
forward, demonstrating the power and flexibility of the frameworks set up in 
|TZ2j and |MaZlj . However, there is an interesting and apparently general 
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difficulty associated with the absence of parabolic smoothing, namely, the 
need to show that nonlinear source terms, defined as the difference between 
nonlinear solutions and solutions of the linearized equations, are quadrati- 
cally small in the X\ norm relative to the X\ norm of initial data, the key 
point being to control the high-derivative norm || • \\jja despite apparent loss 
of derivatives. We discuss this issue in detail in Section 

An analogous issue arises in the nonlinear stability theory in going from 
the strictly parabolic to the hyperbolic-parabolic case, with the difficulty 
again to control an apparent loss of derivatives. This is resolved in |MaZ2[ 
IMaZ3j with an auxiliary nonlinear energy estimate. Here, as there, we ob- 
tain the needed derivative control by auxiliary energy estimates (Prop. 13. 2p . 
However, these are of a rather different type, being local rather than global 
in time, and measuring variational rather than time-asymptotic properties. 
They are also somewhat more delicate, depending strongly on the structure of 
the hyperbolic part of the equations. In particular, they hold in Lagrangian 
but not Eulerian coordinates, indicating the importance of nonlinear trans- 
port effects; see the discussion of Appendix [A] 

A second new feature in the analysis is the incorporation of exponentially 
weighted norms, yielding (I1.12p . We accomplish this by an improved way 
of accounting cancellation, which at the same time greatly simplifies the 
analysis. The key estimate in [TZ2J corresponds heuristically to showing 
that the kernel 

oo 

(1.13) Ky{x,y) :=Y,K y (x,jT;y) 

3=0 

of the ^/-derivative of formal inverse 

oo 

(1.14) (H-e LT )- 1 = E eI * r ' 

3=0 

L := d\ — ad x , converges in L°°(L 2 (x);y) uniformly for < C^ 1 < T < C, 
where 

K(x,t;y) :=cor 1/2 e^-y- a ^/ 4t 

is a convected heat kernel, with a > 0, with decay of JC y (x, y) = JC y (x — y, 0) 
in \x — y\ determining the ultimate decay rate for |U r — U e \. 

Since \\K(x, jT;y)\\ L 2^ = Ct^ 1 ^, convergence in (11.131) cannot be abso- 
lute, but must involve cancellation. This was detected in |TZ2j by using the 
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defining equation K t + aK x — K xx = to write 
(1.15) 



_ r+oo 

^2K y (x,jT;y) ~T~ l K y (x,t;y)dt 

3=1 JT 

r+oo r+oo 

= -a- l T~ l K t (x,t;y)dt + a^T' 1 K w (x,t;y)dt 
Jt Jt 

r+oo 

= a~ 1 T~ 1 K(x,T;y) + a _1 T _1 / K yy (x,t;y)dt, 

Jt 

and observing that K yy (x, t; y)dt is convergent and < C(l + \x — y|) _1 . 

On the other hand, we could just as well have repeated the process to 
estimate 
(1.16) 

r+oo r+oo 

/ Kyy(x,t;y)dt = a~ l T~ l K y (x,T]y) + gT 1 !^ 1 / K yyy (x,t;y)dt, 
Jt Jt 

with | f^°° K yyy (x,t;y)dt\ < C(l + \x — y\)~ 3 ^ 2 - Continuing the process, we 
find that JC y decays to any polynomial order, with asymptotic expansion 

oo 

(1.17) JC y {x, y) = T- 1 a~ j cP y K(x, T; y), 

3=0 

each term of which exhibits exponential decay, but for which we have shown 
only polynomial spatial bounds and not convergence to zero of the trunca- 
tion error of successive finite sums. Moreover, there is still the issue of the 
"continuization error" introduced at the first step of (11.151) . 

In the present analysis, we sidestep these issues, effectively summing to 
all orders expansion (I1.17p . by working at the level of the inverse Laplace 
transform formula 



;i-18) K{x,y;t) = ^e xt K x (x,y)d\ 



used to obtain pointwise bounds for the actual, variable-coefficient problem, 
where K\(x, y) is the resolvent kernel, or Laplace transform of K(x, t; y), and 
T is an appropriate sectorial contour contained strictly in the resolvent set 
of L := d x — ad x and lying in the strictly negative half-plane 3?A < — 770 < 
for |A| > r. 
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Contributions from A > r > are negligible, and likewise from x < y 
(recall, a > 0), and so we may focus on the low- frequency regime |A| < r and 
x > y, where K x (x,y) = c{X)e tl<yX ^ x ^ y ' > with 

(1.19) c(A) = ( Co + Ac 1 + ---) 

/i(A) = A/i(A), fx = —a 1 + di\ + • • • , 

and thus d y K x = Xfi(X)K x . The assumption that T lie in the resolvent set 
implies further that Jfyi < — r] <0forAGrn{|A| < r}. 

Carrying out the sum (11.131) exactly, therefore, within Laplace inversion 
formula (I1.18p . on low frequencies |A < r| and x > y, we obtain a contribution 
to JC y (x, y) of 
(1-20) 

/ y2e XjT d y K x (x,y)d\= I ( V \e x > T ) fi(\)K x (x, y)d\ 

Jrn{|A|<r} j=0 Jrn{\\\<r} v j=0 / 



A 



7rn{|A|<r} v i — e 



x Li(\)K x (x,y)dX 
A 



AT 

rn{|A]<r} V1 — e 



L~i(X)K x (x,y)dX 



-a^T- 1 ! \K x (x,y)\d\ 
irn{|A|<r} 



for some C, r\ > 0, by (11.191) and Jfyt < — ?7 . Not only does this argument 
yield exponential spatial decay bounds, but it turns out to be much eas- 
ier than the previous argument to generalize to the full, variable-coefficient 
problem, thus streamlining and sharpening the analysis at the same time. 

Remark 1.10. Accounting cancellation in this way on the whole solution at 
once is much simpler than breaking the solution into individual components 
and estimating each separately as was done in [TZ2j . Notice further that 
we have essentially performed the inversion (jl,14p directly via the spectral 
resolution formula 

(1.21) (Id - e LT Y x d x = j{l - e^dyGxdX, 
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which has an "elliptic" flavor perhaps somewhat analogous to the spatial 
dynamics point of view. However, in the full problem, there is a boundary 
contribution coming from the J — > oo term in (11.201) for certain (stationary 
eigen-) modes, that is not reflected in formula (11.211) ; see Remark 12.81 



1.5 Discussion and open problems 

In the strictly parabolic case |TZ2] . the analysis of the one-dimensional case 
extended easily to yield a corresponding multidimensional result for shocks 
propagating along a cylinder of finite cross-section with artificial (Neumann 
or periodic) boundary conditions. In the hyperbolic-parabolic case, how- 
ever, our reliance on Lagrangian coordinates in deriving H s energy estimates 
limits us for the moment to one dimension; see Remark IA.1I and the dis- 
cussion of Appendix [A] It would be very interesting (both physically and 
mathematically) to remove this restriction, perhaps by the incorporation of 
nonautonomous effects in the linearized equations as in [D] . 

A further issue when and whether bifurcation condition (T> £ ) is actually 
satisfied in physically interesting situations. It may well be that (T> s ) does 
not occur at all for ideal gas dynamics, nor any other kind of instability. 
However, it seems quite possible that (T> e ) can occur in the richer setting of 
MHD, where Lax-type shock waves are known to be sometimes unstable even 
for an ideal gas equation of state [T] , or for phase-transitional gas dynamics 
with a van der Waals-type equation of state. Numerical investigation of 
(T> £ ) across a range of shock and detonation waves is an important direction 
for future investigation. Numerical evaluation of the Evans function may 
carried out in a well-conditioned fashioned, as described, e.g., in |BrH [Br2, 
iBrZl IBDUI IHuEl IBHR] . 

The analysis of the present paper serves also as a stepping-stone to the 
closely related but more complicated detonation case. With suitable elabo- 
ration, the argument extends to that case [TZ3j . yielding Hopf bifurcation, 
or "galloping instability" of viscous strong detonation waves in one dimen- 
sion. As in the shock case, this result of course is subject to verification of 
the spectral bifurcation hypothesis (T> e ). For detonation solutions of react- 
ing gas dynamics, galloping is expected to occur, and Evans function results 
|LyZl| |LyZ2| show that instability if it occurs must be of a (possibly degen- 



erate) Hopf type spectral configuration. 

Finally, we mention the problem of determining stability of the bifurcat- 
ing time-periodic waves whose existence is established in Theorem II .81 In the 
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absence of a spectral gap, our method of analysis does not directly yield sta- 
bility as in the finite-dimensional ODE case, but at best partial information 
on the location of point spectrum associated with oscillatory modes, with 
stability presumably corresponding to the standard condition de/dr > 0; see 
|SS] for a corresponding result in the strictly parabolic case. The hope is 
that one could combine such information with an analysis like that carried 
out for stationary waves in [ZHj, IMaZlj , adapted from the autonomous to the 
time-periodic setting: that is, a generalized Floquet analysis in the PDE set- 
ting and in the absence of a spectral gap. This seems a particularly exciting 
direction for further development of the theory. 

Plan of the paper. We begin in Section [2] by recalling the linearized 
bounds of |MaZl[ IZ2| ITZ2j . and establishing the new bound described in 
f ll.20p . In Section [31 we carry out the nonlinear energy estimates controlling 
higher derivatives in nonlinear source terms. In Section HI we prove the main 
theorem, carrying out the bifurcation analysis following the framework of 
|TZ2] . Finally, in Section we show that our results apply to gas- and 
magnetohydrodynamics. 

Acknowledgement. Thanks to Bjorn Sandstede and Arnd Scheel for 
their interest in this work and for stimulating discussions on spatial dynam- 
ics and bifurcation in the absence of a spectral gap. B.T. thanks Indiana 
University for their hospitality during the collaborative visit in which the 
analysis was carried out. 

2 Linearized estimates 

Assuming (Al)-(A2), (H0)-(H4) (alternatively, the weakened hypotheses of 
Remark I3.5j) . let L £ as in the introduction be the linearized operators (II. 3p . 
and \±{e) the crossing eigenvalues of (T> £ ). 

Lemma 2.1 (TZ2). Associated with eigenvalues \±{e) of L e are right and 
left eigenfunctions <fi e ± and <p e ± G C k (x, e), k > 2, exponentially decaying in 
up to q derivatives as x — > ±oo, and L e -invariant projection 



onto the total (oscillatory) eigenspace S £ := Span{0±} ; bounded from L q or 



(2.1) 




j=± 
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B 2 to W k ' p fl X 2 for any 1 < q,p < oo. Moreover, 
(2.2) £ ± = d x & ± , 

with $ £ G C k+1 exponentially decaying in up to k + 1 derivatives as x — > ±00. 

Proof. From simplicity of A±, and the fact |MaZlj that they are bounded 
away from the essential spectrum of L E , we obtain either by standard spectral 
perturbation theory |Katj or by direct Evans-function calculations [GJltlGJ2| 
IMaZHlZ2] that there exist A±(-), <j>±(-) G L 2 with the same smoothness C (e), 
k > 5, assumed on F. The exponential decay properties in x then follow by 
standard asymptotic ODE theory; see, e.g., \GZ\ IZ2l IZ3j . Finally, recall 
the observation of [ZHJ that, by divergence form of L e , we may integrate 
L £ (j) = from x = —00 to x — +00 to obtain A f^°° (f){x)dx = 0, and 
thereby (since A ± 7^ by assumption) 

/+00 
<p±{x)dx = 0, 
-00 

from which we obtain by integration (12.21) with the stated properties of $±. 
From f)2.2p and representation (12. ip . we obtain by Holder's inequality the 
stated bounds on projection II. □ 

Defining IF := Id — IF, S e := RangelF, and L £ := L £ U £ , denote by 

(2.4) G(x,t;y) := e L ^6 y (x) 

the Green kernel associated with the linearized solution operator e Lt of the 
linearized evolution equations u t = L e u, and 

(2.5) G(x,t;y) :=e iet US y (x) 

the Green kernel associated with the transverse linearized solution operator 
e Le *II. By direct computation, G = O + G, where 

(2.6) 0(x,t;y) := e { ^ £)+iT{£))t <j> + (x)ft + {y) + e (7(£) - tr(e)) V-(x)0*_(i/). 
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2.1 Short time estimates 

Lemma 2.2. For < t < T , 1 < p < oo, rj > 0, and some C = C(T), 



(2-7) ||e Le 7lUp, ||e^n/|| LP <q|/|| 



LP- 



(2.8) ||e"<*> / e L ^d x f\\ LP , \\e^ / e L ^fld x f\\ LP < C\\e^f\\ 



\LP, 



[2.9) We"® J d e , t e L °%f\\ LP , \\e^ J d £ , t e L ^Ild x f\\ LP < C\\e^f\\ 



W 2P : 



where (x) := (1 + l^ 2 ) 1 / 2 . 

Proof. From standard C° semigroup bound \e Let \i jP ^i jP < C and properties 
e L e t = e L E t u + e i s tfi and ||n/|| LP < \f\ LP} we obtain (j277l) . Likewise, we 
may obtain integrated bounds (12.81) (r = 0) using the divergence form of L £ , 
by integrating the linearized equations with respect to x to obtain linearized 
equations U t = C £ U for integrated variable 

/X 
u(z,t,e)dz, u(-,t) :— e Lst d x f, 
-oo 

with linearized operator C £ := — A £ (x)d x + d x of the same parabolic form 
as L e , then applying standard C° semigroup estimates (alternatively, more 
detailed pointwise bounds as in |MaZlj ) to bound 

\\e^U(-,t,e)\\ LP = \\e^e c ^f\\ LP < C\\e^ f\\ LP , 

the e L *Il bound then following by relation e Le * = e Le< n + e ie *Il together with 
|| e ij<*> J Ud x f\\ LP < \ f\ LP . We obtain (JHED (r > 0) by the change of variables 
V = Ue v ( x \ which, since a := satisfies \(d/dx) k a\ < C\a\ for k > 0, 
converts the linearized equations to an equation with the same principal part 
plus lower-order terms with bounded coeffients. 

Finally, (I2.9P follows from dte 1 ^ 1 = e L ^- ) *L(e) and the variational equa- 
tion (d t - L)d £ U = (d £ L)U, U(t) := e L ( £t U , together with \\L(e)U \\ LP , 
\\L(e)U\\ LP <C\\U \\ W 2, P . □ 
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2.2 Pointwise Green function bounds 



We now develop the key cancellation estimates analogous to (jl.20p of the 
introduction, adapting the pointwise semigroup methods of [ZH] IMaZl[ IZ2] 
to the present case. Our starting point is the inverse Laplace transform 
representation 

I prj-i-ioo 

(2.10) G(x,t;y) = — P.V. / e xt G x (x,y)d\, 

7] > sufficiently large, established in |MaZlj . 

Deforming the contour using analyticity of G\ |MaZl] across oscillatory 
eigenvalues A±(e) we obtain G = G + O, where O, defined in (|2.6|) . is the sum 
of the residues of the integrand at A±, and, for u, r > sufficiently small, 

G{x,t;y) = ^-£e xt G x (x,y)d\ 

p—u—ri p—u+ioo 

P.V.(/ +/ )e xt G x (x,y)d\ 

J —u—ioo J —u+ri 



(2.11) 1 

=: & + G 11 , 



where 7 is the counterclockwise arc of circle dB(0, r) connecting —v — ri and 
—u + ri, and G 1 as in |Z3J is the low-frequency and G 11 the high-frequency 
component of G. Define associated solution operators S ,/ (t,e) and S n (t,e) 
by 

p+00 

(2.12) SP(t)f(x):= / &{x,t;y)f{y)dy 



and 

S : = S 1 + S 11 . 

Supressing the parameter e, denote by af, r^, W the eigenvalues and 
associated right and left eigenvectors of A £ ± = F u {u e ±1 e). Following [MaZlj . 
let a*(x), j = 1, . . . , J < (n — r) denote the eigenvalues of 

A* := An — 5 2 i-B 22 1 ^4i2 = A n . 

Let L(X, Y) denote the space of bounded linear operators from Banach space 
X to Y, equipped with the usual operator norm | • \l(x,y)- 
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Proposition 2.3. Under assumptions (A1)-(A2), (H0)-(HJ h ), (D £ ) (alter- 
natively, those of Remark \3 . 5\) , for < q, r, s, < 2q + s < A, and v > 0, 
(2.13) 

j 

d?d r t d s x G n (x,t;y) = (E°( e i^("!/) + ^ (N+t) ))^ 

p<(2q+r)+s j=l 

Proof. The case (q = r = 0) is immediate from the bounds of [MaZlj . Deriva- 
tives with respect to e may be converted using the variational equation 

(2.14) (\-L)(d £ G x ) = (d £ L)G x 

into two spatial derivatives, thus extending to case (r = 0). Likewise, time- 
derivatives, appearing as powers of A at the level of the resolvent kernel, have 
the effect of single spatial derivatives on "hyperbolic" blocks, and double 
spatial derivatives on "parabolic" blocks, in the notation of [MaZlj , yielding 
by the same estimates as in case (q — r — 0) the asserted bounds. This 
completes the proof of the general case. □ 

Corollary 2.4. Under assumptions (A1)-(A2), (HO)-(H^.), (D e ) (alterna- 
tively, the weakened assumptions of Remark \ 3. 50 , YlJLo S 11 (jT) converges 
uniformly and absolutely in operator norm | • \uxi,Xi), fore sufficiently small 
and T in any compact set, for any < 77 < v, to a limit that is C 1 in (e, T) 
with respect to the operator norm L(Bi, Bi). 

Proof. Straightforward computation using (12.131) . □ 

Proposition 2.5. Under assumptions (A1)-(A2), (H0)-(H4), (D e ) (alter- 
natively, the weakened assumptions of Remark \3.5\) . YlJLo^ijT) converges 
uniformly (but conditionally) in the operator norm | • \l(x 2 ,b 1 ), for e suffi- 
ciently small and T in any compact set bounded away from the origin, with 
limiting kernel satisfying for < q, r, s, < 2q + s < 4, v > 0, 

00 

(2.15) &&ar m Gl(x,jT; y) = 0(e'^ + e -"M). 

3=0 

Proof. (Convergence). Since (by Sobolev embedding) d x \\e 2ri ^ x ' ) ■ C X 2 , 
it is sufficient to show convergence in L(d x \\e 2r1 ^ ■ ||ioo,L 2 ), or, equivalently, 
convergence in L(||e 2,? ^ ■ ||loo,L 2 ) of the operator YlJLo^ 1 ^ with kernel 

j 

^2Gl(x,jT;y), 
3=0 
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where, by fl2TTTj) . 

(2.16) Gl(x,jT;y) = -L I e ^ T d y G x (x,y)d\, 

hence (summing under the integral as described in the introduction) 



(2.17) 




=: / + II. 



Taking y < for definiteness, recall from [MaZl| IZ2j that, for A suffi- 
ciently small, G\(x,y) can in the Lax or overcompressive case be expanded 
analytically at A = as a sum of "excited" terms 

(2.18) A-V m (*)^V^, A* = -(aj) -1 + 0(A), 

where 0j(x) = 0(e _i/ ' x '') are stationary modes of the linearized equations; 
"scattering terms" 

(2.19) r p-T e X^(X)y + Xflf(X)^ ~± = _ (a ± } -l + 0(A) . 

and error terms 

(2.20) o(e' u{]x]+lv]) ) 
and 

(2.21) o(\e x ^ (x)y+x ^f Wx ), 

with symmetric expansions for y > 0. These bounds may be converted by 
the Riemann Saddlepoint (~ Stationary Phase) estimates to the pointwise 
Green function bounds of [ZHl IMaZIL E2\ . 
Taking y-derivatives, we find that 

( T ^ T )d y G x = -T-^X- 1 + 0(l))d y G x 
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expands as (Ta^T ) _1 times the same excited and scattering terms (12.181) and 
( I2.19P , plus error terms of the same order (I2.20p , (I2.2ip , plus new pole terms 
of order 

(2.22) \- 1 0(e- v{ ^ + ^). 

Thus, by the same Riemann saddlepoint estimates used to bound &{x, t; y) 
in [ZH[ IMaZl[ IZ2] , we find that the contribution to term II of all except 
the new terms (12.221) satisfies exactly the same bounds as Gr(x, (J + 1)T; y), 
expanding for t = ( J + 1)T as the sum of excited terms 



(2,3, (r.nwhfe|)-^(^)), 

(3^ > 0, scattering terms 

n 

5^(1 + t)-i/2 (e-(-^ a fe "*) 2 / Mi e-^ + ) 
fc=i 

(2.24) + 2^ ^{l«fc*l>bl}^ + ^ e e 

a~>0, a~<0 

+ E X { | a - t |> b | } (l + t)- 1/2 0(e-(— ^(H»/-*-|)m* e -^ 



bounded by convected heat kernels, where x + (resp. x ) denotes the positive 
(resp. negative) part of x, and a neglible error 

(2.25) 0(e- u{lx - yl+t) ). 

Terms (12.221) by the same argument give a time-independent contribution of 

(2.26) ^fm(x)9m(y) = 0(e-^ + W) 

up to negligible error (I2.25p . 

Noting that all terms except (12.231) and (I2.26P decay in L 2 (x) at least at 
Gaussian rate (1 + 1)^ 1 ^ 4 , independent of y, while kernels (12.231) integrated 
against an exponentially localized function converges exponentially to 

= 0( e -"M), 
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we find that II converges as J — > oo to a kernel with the claimed bound 
Q^ e -v\x\ _|_ e v\x-y\y Term I (independent of J) by the elementary bound 
\d y G x (x,y)\ < Ce- u \' x -y\ for A bounded away from essential spectrum of L, 
and the fact that 7 is bounded away from both A = and the essential 
spectrum of L, satisfies trivially the bound / = 0(e~ v ^ x ~ y ^), completing the 
proof of convergence while at the same time establishing (12.151) for q,r,s — 0. 

(Derivative bounds). The remaining bounds (I2.15p (g. r. s ^ 0) follow 
easily, by the observation that the limiting summands in term II are inde- 
pendent of T and satisfy the same bounds after x- or e- differentiation as 
before, while J, by (either using the variational equations (12.141) to convert 
e-derivatives to two spatial derivatives as above, or estimating directly as in 
[TZ2] ) 

d?d s x d y G x (x,y) = 0(e-^-y\) 

and 

i*(t^)| £C 

for A G 7 likewise satisfies the same bounds after as before. □ 

Remark 2.6. Note that, except for the excited term (12.231) . the estimates on 
II are rather crude compared to the detail in [MaZl, Z2J, while the estimate 
for / is trivial. This represents a considerable simplification from |TZ2j . 

Remark 2.7. In the undercompressive case, there appear in the expansion 
of G\ additional terms of order X^e'^^ and e A ^ ±(A)x e"^ |y| that do not 
fit in our analysis. As discussed in [ZHl IZ2j . this case is essentially different. 

Remark 2.8. Note that there is indeed a contribution from the limiting 
term // at infinity, as mentioned in Remark II. 101 

Corollary 2.9. Under assumptions (A1)-(A2), (H0)-(H4), (V e ) (alterna- 
tively, the weakened assumptions of Remark \3.5\) . Y^jLo^tiT) converges uni- 
formly in operator norm | • \l(x 2 ,b 1 ), for e sufficiently small and T in any 
compact set bounded away from the origin, to a limit that is bounded in 
I • |i(x 2 ,Xi)ni(B 2 ,-Bi) an d C l (e,T) in \ ■ \l{b 2 ,Bi) on the subspace X 2 . 

Proof. By Corollary 12.41 and Proposition 12. 5[ it is sufficient to establish the 
second statement for S 1 using bounds (I2.15p . Further, since d x L x C B 2 and 
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(by Sobolev embedding) d x \\e 2r) ^ • ||ioo C X 2 , it is sufficient to show bounded- 
ness in L(d x L 1 , L 2 ) and L(d x \\e 2r! ^ ■ ||i«>,Xi), or, equivalently, boundedness 
in -^(L 1 , L 2 ) of the operator 9(£, T) 9 ^°1 S I d x , with kernel 

oo 
j=0 

q = 0, 1, and boundedness in L(||e 2r? ^ ■ ||x«>, ||e' ? ' a: ' ■ ||£,«,) of the operators 
3£ Sjlo ^ 7 ^' < r < 3 with kernels 

3=0 

both routine consequences of bound (12.151) . □ 
From Corollary 12.91 we obtain the following important conclusion. 

Proposition 2.10. Under assumptions (A1)-(A2), (H0)-(H4), (D e ) (alter- 
natively, the weakened assumptions of Remark \3. 5\) . fore sufficiently small, 
and < Tq < T < T%, (Id — S(s, T)) has a right inverse 

(Id-S(£,T))- 1 GL(X 2 ,X 1 ) 

that, restricted to its domain X 2 , is C 1 in e, T in the L(B2,Bi) norm. 

Proof. The first assertion follows by a standard telescoping sum argument, 
setting (Id — S^T)) -1 := Yl'jLo ^( £ > -O" 7 ' the second by the final assertion 
of Corollary 12. 91 See [TZ2j . Section 2, for further details. □ 

3 Nonlinear energy estimates 

We next carry out .£P-energy estimates on the perturbation equations 
(3.1) U t -L £ U = Qe(U, U x ) x 

of (11.21) about U £ , to be used to control higher derivatives in the fixed point 
iteration used to carry out our Lyapunov-Schmidt reduction and bifurcation 
analyses. 
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By standard energy estimates, we have H s well-posedness, 

\\U(;T)\\ H s<C\\U \\ H s 

for T uniformly bounded. Likewise, we have a formally quadratic linearized 
truncation error \Q\ = 0(\U\(\U\ + \U X \) for \U\ < C. Our goal in this 
section, and what is far from evident in the absence of parabolic smoothing, 
is to establish a quadratic bound on the linearization error. 

(3.2) \\U(;T)-e L ° T U \\ H s <C\\U \\ 2 H .. 

The corresponding bound does not hold for quasilinear hyperbolic equations, 
nor as discussed in Appendix |A] for systems of general hyperbolic-parabolic 
type, due to loss of derivatives. However, it follows easily for systems satis- 
fying assumptions (Al)-(A2). 



3.1 H s linearization error 

By the strong block structure assumption (Al), we may write (13. ip more 
precisely as 



(3.3) U t + {A £ U) X - {B £ U X ) X = d x 





qf{U){U,U) + ql(U)(U,d x U 2 ] 



where B £ = ^ ^.J and qj(-) are C k 2 bilinear forms. Appealing to (A2), 



we may differentiate i times and multiply by 
to obtain 

(3.4) A ^diU t + A%U X + -(B £ d e x +1 U) x = G t + A°> £ d x 



A?f 
A% 







where A°' £ is symmetric positive definite, Af 1 is symmetric, B £ = ^ ~A 
with b £ symmetric positive definite, and 

||e,|| L2 = o(||t/|M, 

qt = di{ql{U){U,U)+ql{U){U,d x U 2 )). 
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Lemma 3.1. For t > 1, \\qA\ H t < C(\\U\\ H e + \\U 2 \\ H t + x) (\\U\\ H t + ||C/||^). 
Proof. Standard application of Moser's inequality; see, e.g., \Ta\ IZ3j. □ 

Proposition 3.2. Assuming (Al), (A2), (HO), for 1 < s < k - 1, < T < 

T uniformly bounded, some C = C(T ) > ; and U satisfying (13. ip with 
initial data U(-,0) = U sufficiently small in H s , 

(3-6) \M;T)\\ 2 HS + [ T \\U 2 (;t)\\ 2 Hs+1 dt<C\\U \\ 2 HS , 

Jo 



(3.7) \\U(;T)-e LT Uo\\Hs<C\\U \\ 2 H s. 
Proof. By symmetric positive definiteness of A 0,£ , 

S(U):= (1/2) J2(diU,A^diU) 
i=i 

defines a norm equivalent to || • i.e., S(-) 1 ^ 2 ~ || • \\h s - Applying (13. 4p . 
we find that 

d t £(U) = -(d s x U,A s d s x +1 U) + (d s x U,d x (B*d x U)) 
+ (d s x U 2 ,A^d x q s )+0(\\U\\ 2 H3 ) 
(3 g) = (d s x U, [l/2)AldtU) + (d s x +1 U 2 , tfd s x +1 U 2 ) 

- (d s x +1 U 2 ,A^q s ) - (d s x U 2 ,d x A° 22 £ q s ) + 0(\\U\\Hs) 
<-9\\U 2 \\ 2 hs+1 +0(\\U\\hs) 
+ 0(\\U 2 \\ Hs+ i) (||E7|| H , + \\U 2 \\ H s + i)(\\U\\ H s + \\U\\ S HS ) 

for some 6 > 0. So long as remains sufficiently small, this gives 



dtS{U) < -{6/2)\\U 2 \\ 2 HS+1 + 0{\\Uf HS 
<-(e/2)\\U 2 \\ 2 Hs+1 +C£, 



(3-9) . „,MM,„„2 



from which (13.61) follows by Gronwall's inequality, in the form 

S(U(T)) + (9/2) f \\U 2 \\ 2 Hs+1 (t)dt < C 2 S(U ). 
Jo 
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To obtain (13.71) . observe that V(-,t) := U(-,t) — e Lst U satisfies 

(3.10) E t + (A £ E) X - (B^E X ) X = d x ( gf(f/)(f/? v) + ° q e {um dxU2) ) , 
with E(-, 0) = 0, and thus 

(3.11) A°^d e x E t + A"d e x E x + -(B% +1 E) X = Of + A°^d x , 

similarly as in (13.41) . with \\Qf \\ L 2 = 0(\\E\\ H i). Thus, calculating as in (13. 8p . 

we obtain 

(3.12) 

d t £(E) = -(d s x E,A £ d: +1 E) + (d x E,d x (&d s x E)} 

+ (d s x E 2 ,A°^d x q s ) + 0(\\E\\ H s(\\E\\ H s + \\U\\%.) 

= (d x E, (l/2)Ald°E) + (d s x +1 E 2 ,b*d s x +1 E 2 ) 

- (d s x +1 E 2 ,A°^q s ) - (d s x U 2 ,d x A%%) +0(\\E\\ H .(\\E\\ H . + \\U\\ 2 Hs ) 

< -9\\E 2 \\ 2 HS+1 + +0(\\E\\ H s(\\E\\ H s + \\U\\%.) 
+ 0(\\E 2 \\ H s +1 )(\\U\\ H s + ||[/ 2 ||^+i)(||C/||^ + \\U\\ S HS ), 

which, so long as H^H/p remains sufficiently small, gives 

d t £(E) < -9\\E 2 \\ 2 HS+1 + C\\Ef HS + C\\U\\ 2 H . (\\U\\ 2 HS + \\U 2 \\ 2 Hs+1 ) 

and thus, using < C , ||?7o||h s by (13.61) . 

d t S(E) < +CS(E) + C 2 \\U \\ 2 H .(\\U\\ 2 H . + \\U 2 \\ 2 HS+1 ). 
Applying Gronwall's inequality, and using -E'(O) = 0, we thus obtain 

£{E{T)) < C 3 \\U \\ 2 H s [ (\\U\\ 2 HS + \\U 2 \\ 2 HS+1 )(t)dt, 

Jo 

yielding the result by and £(E) ~ \\E\\ 2 Hs . □ 

Remark 3.3. Note that the above, finite-time estimate is considerably sim- 
pler than the "Kawashima-type" global-in-time estimates used in the nonlin- 
ear stability analysis |MaZ3j |Z2| |Z3] , which require also additional assump- 
tions of, among other things, symmetrizability and "genuine noncoupling" of 
A% and B e ± . 
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3.2 Weighted variational bounds 

Define now 

(3.13) N(Uo, e, T) := U(-, T) — e^ T U . 

where U satisfies ( 13. D with initial data Uq. 



Proposition 3.4. Assuming (Al), (A2), (HO), for 1 < s < k - 1, < T < 

T uniformly bounded, some C = C(T ) > 0, Af is uniformly bounded in X 2 
and C l in B 2 with respect to Uq, T, e for \\Uq\\x 1 , \e\ sufficiently small, with 

||^o,e,r)||x a <C7||^o||^, 
(3.14) \\d Uo N(Uo,e,T)\\ X2 <C\\U \\ Xl , 

\\d T;E M(Uo,e,T)\\ B2 <C\\U \\ Xl . 

Proof. Similarly as in the proof of fl2.8j) (t > 0), we may make the change 
of variables V = Ua(x), a(x) := e n ^ , (x) := (1 + |x| 2 ) 1//2 , to convert both 
(13.11) and the corresponding linearized equations to equations for which the 
same principal part plus lower-order terms with bounded coeffients and the 
nonlinear part is a^ 1 times the same principal part plus bounded factors 
times lower-order terms, thus recovering \\U(-,T) — e LT Uo\\x 1 < C||£7o||xi> or 



(3.15) \\tf(U ,e,T)\\ Xl <C\\U Q \ 



by the same argument used to obtain (I3.6p . Here, we are using the fact that 
a satisfies both < C\ and \(d/dx) k a\ < C\a\ for k > 0. Bound (I3.14p (i) 
then follows by (13.151) and the Duhamel formulation 



r 



Af = J v z m{T ~ t)d *)Q( u M{t) d t, 

bounds (12.81) with p = oo and r = 2, and Sobolev embedding 

\\e^U(;t)\\ L ~, \\e^d x U(;t)\\^ < \\e^U(-,t)\\w^ 

<C\\e«MU(;t)\\H>. 

Bound (I3.14p (ii) follows, similarly, by repeating in weighted norm Xi the 
arguments for (13.61) and (13 .7p with U replaced by W := U 1 — U 2 , where f/ J 
are solutions of (13.11) with different initial data Uq and Uq. Bound (I3.14p (iii) 
follows by (13.141) (i) together with the equation for F := W — e L ^ T W Q . □ 
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Remark 3.5. It is readily checked that bounds fl3.6p . fl3.7p . f)3.14p hold also 
under the weakened assumptions described in Remark 11.31 m the favorable 
variables (U\, W 2 ), by rearranging (11.21) as 



(3.17) 



U ± \ + (A n {U) *\ [U x 



W 2 J \ * * J \W 2 



\(u x 

(dW 2 /dU 2 )P(U) \W 2 J^\0(\U x ^ 



similarly as in [Z2| IZ31 IGMWZj and carrying out energy estimates as be- 
fore. These then imply corresponding estimates in the original variables 
(Ui,U 2 ), by the observation that, for U satisfying the linearized equations 
in (Ui,U 2 ) coordinates, (Ui, (dW 2 /dU)(U)U) satisfies the linearized equa- 
tions in (U\, W 2 ) coordinates, and, for U = U — U satisfying the nonlinear 
perturbation equation in (U\,U 2 ) coordinates, (Ui,W 2 (U)) — {Ux,W 2 {U)) = 
(Ui, (dW 2 / ' dU)(U)U) +0(|?7| 2 ) satisfies the nonlinear perturbation equation 
in (Ui,W 2 ) coordinates. 

4 Bifurcation analysis 

We now carry out the bifurcation analysis following the framework of [TZ2J, 
with a few slight modifications to simplify the analysis. 

4.1 Construction of the period map 

Given a solution U of (11.21) . define the perturbation variable 

(4.1) U(x, t, e) := U £ (x, t) - U £ (x), 
satisfying nonlinear perturbation equations 

(4.2) U t - L e U = Q E (U, U x ) x , U(x, 0, e) = U (x, e). 

By dan, <m, 

(4.3) Q e (U,e)=0(\U\ 2 + \U\\U x \) 

so long as U satisfies a uniform L°° bound: in particular, for ||Z7||xi ^ C. 



26 



Decomposing 
(4.4) U = w l (j) + + w 2 4>- + v, 

where Wi4> + + W2(p~ G £, v :— UU G S, and coordinatizing as (w,v), we 
obtain after a brief calculation 



if 

(4.5) 



2% $j) «+*.<«,«,«), 



u = L £ t> + N v (w,v,e), 

where L = LU and 

, , N w>1 <f> + + jV w>2 0_ = IIQ e (tf, U X ) X1 

(4.6) 

N v = UQ e (U,U x ) x . 

Lemma 4.1. N w is C k , k > 2, from (w,v) G R 2 x L q to R 2 , for any 

< q < oo, wzt/i 

(4.7) |i\g<C(M 2 + IMlU 

iV„ C fc ; jfe > 2, /rom (w, u) G R 2 x X l to X' 2 := d x {f : \\e 2 ^ f\\ H i < +00} 
and, for (\w\ + IMIxJ < C, C 1 /rom G R 2 x Si to dx^L 1 ), with 

(4.8) ||iV„||^ <C(|^| 2 +\\v\\ 2 Xi ) and pA^k^) < C(M + \\v\\ Xl ). 
Proof. Direct calculation, using H4.3|) and the H-bounds of Lemma 12.11 □ 

Proposition 4.2. For < t < T , any fixed Ci, T > 0, some C > 0, and \a\, 
||&||xi; kl sufficiently small, system ( 14. 5p with initial data (wq,Vo) = (a,b) 
sufficiently small in R 2 x X% possesses a solution 

(4.9) {w,v){a,b,e,t) G R 2 x X 1 

that is C k+l in t and C k in (a,b,e), k > 2, with respect to the weaker norm 
B\, with 

f4im - C\bf Xx < \w(t)\ < C(\a\ + \\b\\ 2 Xl ), 

1 ' 1 IK*)lk< Wlk + M 2 ), 
and 

(4.11) \D M (w,v)(t)\ R2xBl ^ Bl <C. 
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In particular, for \\b\\xi < CiM, all < t < T , 

(4.12) \\v(t)\\ Xl < C\w(t)\. 
Likewise, for \\b\\x 1 < Ci |a| 2 , all <t <T, 

(4.13) \\v(t)\\ Xl <C\w(t)\ 2 . 

Proof. Existence and uniquess follow by a standard Contraction-mapping ar- 
gument, using a priori bounds (13.141) . which also imply (14.101) (by decoupling 
of linear parts) and (14. lip . Combining (I4.10l) (i)-(ii). we obtain evidently 
ff47T2|) and fl4TT3|) for \a\ sufficiently small. □ 

Setting t = T in (14. 9ft and applying Duhamel's principle/variation of 
constants, we may express the period map 

(4.14) (a, b, e) -> (a, b) := (w, v)(a, b, e, T) 
as a discrete dynamical system 

a = R(e, T)a + Ni(a, b, e, T), 

(4.15) 

b = S(e,T)b + N 2 (a,b,e,T) 

with e, T e R 1 , a, Ni E M 2 and b e B 1 , N 2 e B 2 , where 

R(e,T) := e 7(£)W ' £) , 
(4-16) . ' 

S(e,T) := e L * T ^ 

are the linearized solution operators in w, v and 

r T 

#i(a,&,e,r) := / e^ e)(T ' s) N w (wv, e)(s)ds, 
(4.17) Jo 

N 2 (a,b,e,T) := / e L ^ T ~ s) N v {w, v, e)(s)ds 



the differences between nonlinear and linear solution operators: equivalently, 

Ni (a, b, e, T) = HA/"(ai0+ + a 2 0_ + 6, e, T) 
(4.18) ~ , 

N 2 (a, 6, e, T) = (Id - U)Af(a l( j) + + a 2 <\>„ + 6, e, T), 
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where M is the linearization error M defined in ( I3.13P . 

Evidently, periodic solutions of (14. 5 p with period T correspond to fixed 
points of the period map (equilibria of (I4.15P ) or, equivalently, zeroes of the 
displacement map 

A^a, b, e, T) := (R(e, T) - Id) a + N^a, 6, e, T), 
l ' J A 2 (a,6,e,T) := (S(e,T) -Id)b + N 2 (a,b,e,T). 

4.2 Lyapunov— Schmidt reduction 

We now carry out a nonstandard Lyapunov-Schmidt reduction following the 
"inverse temporal dynamics" framework of [TZ2J, tailored for the situation 
that (S(e, T) — Id) is not uniformly invertible, or, equivalently, cr(L) is not 
bounded away from {jn/T}, j G Z. In the present situation, L has both an 
^-dimensional kernel (Lemma 14.51 below) and essential spectra accumulating 
at A = 0, and no other purely imaginary spectra, so that (S(e,T) — Id) = 
(e LT II — Id) inherits the same properties; see [TZ2] for further discusssion. 

Our goal, and the central point of the analysis, is to solve A 2 (a, b, e, T) = 
for b as a function of (a, e, T), eliminating the transverse variable and reducing 
to a standard planar bifurcation problem in the oscillatory variable a. A 
"forward" temporal dynamics technique would be to rewrite A 2 = as a 
fixed point map 

(4.20) b = S(e,T)b + N 2 (a,b,e,T), 

then to substitute for T an arbitrarily large integer multiple jT. In the 
strictly stable case ^ka(L) < —r] < 0, \S(e, jT^nx-Ljd) < 1/2 for j suffi- 
ciently large. Noting that N 2 is quadratic in its dependency, we would have 
therefore contractivity of (I4.20p with respect to b, yielding the desired re- 
duction. However, in the absence of a spectral gap between cr(L) and the 
imaginary axis, \S(e, jT)\ L ( Xl ,Xi) does not decay, and may be always greater 
than unity; thus, this naive approach does not succeed. 

The key idea in [TZ2J is to rewrite A 2 = instead in "backward" form 

(4.21) b = (Id - S(e, T))~ 1 Ar 2 (a, b, e, T), 

then show that (Id — S)^ 1 is well-defined and bounded on RangeiV 2 , thus 
obtaining contractivity by quadratic dependence of iV 2 . Since (Id — S)^ 1 N 2 
is formally given by Yl'jLo S^N 2 this amounts to establishing convergence: a 
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stability /cancellation estimate. Quite similar estimates appear in the non- 
linear stability theory, where the interaction of linearized evolution S and 
nonlinear source N 2 are likewise crucial for decay. The formulation (14.211) 
can be viewed also as a "by-hand" version of the usual proof of the standard 
Implicit Function Theorem |TZ2j . 

Lemma 4.3. Assuming (Al), (A2), (H0)-(H4), (D e ), Ni is quadratic order 
and C l from M 2 x L q x M 2 — > IR 1 for any 1 < q < 00, and N 2 is quadratic 
order from M 2 x X x x M 2 -> X 2 and C 1 from tfx^xl 2 ^ B 2 for \\b\\ Xl 
uniformly bounded, with 

f4 22) \Ni(a, b, e,T)\, [d^N^a, b, e, T)\ L ^ m < C{\a\ + ||6|| Bl ) 2 , 

{daN^b^^lc^^ + IM^M.e^k^) < C(\a\ + \\b\\ Bl ), 

(4.23) 

||JV 2 (a, 6, e, T) \\ X2 , \d E , T N 2 (a, b, e, T) \\ m *,B 2 ) < C(|a| + ||6|| Xl ) 2 , 
\d a N 2 (a,b,e,T)\ c{R2yB2) + \d b N 2 (a, b, e, T)| £(Sli%) < C(|a| + \\b\\ Xl ). 

Proof. Bounds (I4.22p follow from representation (I4.18p . variational bounds 
(I3.14p of Proposition I3.4[ and the Il-bounds of Lemma 12.11 from which we 
likewise obtain 

\\N 2 {a,b,e,T)\\ Xl , \d E , T N 2 {a,b,e,T)\\ m 2 iBl) < C{\a\ + \\b\\ Xl ) 2 , 
\d a N 2 (a, b, e, T)\ m 2 tBl) + \d b N 2 (a, b, e, T)\ £{BljBl) < C{\a\ + \\b\\ Xl ). 

The remaining bounds 

\\N 2 (a, b,£,T)\\ dxe -2 V (x) H i, \d ejT N 2 {a,b,e,T)\\ L ^2 j9a:e -2f,{ X ) H ^ < C(\a\ + \\b\\ Xl 

\d a N 2 (a,b,e,T)\ C (M?,d x Li) + \d b N 2 (a, b, e : T)\ C {B^,d^) < C(\a\ + \\b\\ Xl 

follow easily from Duhamel representation (I4.17p (ii). bounds (I4.10l) (ii). (14.111) . 
and (I4.8p . and bounds (12. 8p . (12.91) on the linearized solution operator e ft. 

□ 

Corollary 4.4 ( [TZ2] ). Assuming (Al), (A2), (H0)-(H4), (V £ ), 

(4.24) A 2 (a, b, e, T) = (a, 6, e, T) G M 2 x Xi x ]R 2 , 
equivalent to 

(4.25) 6 = (Id - 5(0, 6, e, T)) _1 JV 2 (a, b, e,T) +lu 
for 

(4.26) G Ker(Id — 5*(a, 6, e, T)) PI X x . 
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Proof. Applying to the left of (14.241) the right inverse (Id — S(e, a,b)) 1 given 
by Proposition 12.101 we obtain by Lemma 14.31 

b := (Id - S(e, a, 6)) _1 (Id - S(e, a, b))b 
= (Id - S(e, a, 6))- 1 A^ 2 (e, a, b) e X v 

Observing that b — b belongs to Ker(Id — S(e, a,b))n Xi by the right inverse 
property, we obtain (I4.25p . Conversely, (14.251) implies (I4.24p by application 
on the left of Id — S(a, b, e, T). □ 

Lemma 4.5. The kernel of (Id — S(e, T)) is of fixed dimension £ as in (H4), 
is independent of T, and has a smooth basis to = (c^i . . . uij (e) . In the 
Lax case, it is generated entirely by translation invariance. 

Proof. This follows by the corresponding properties of KerL assumed in (H4). 

□ 

Corollary 4.6. Assuming (Al), (A2), (H0)-(H4), (V £ ), the map 

(4.27) T(a, 6, e, T, a) := (Id - S(e, T))- 1 A^ 2 (a, 6, e, T) + w(e)a, 

(Id — S)^ 1 : X% — > X2 as defined in Proposition \2.1Ck is bounded from R 2 x 
Xi x R 2+i -> X x and C 1 from R 2 x B x x M 2+£ -> Xi /or |a| bounded and 
\ a \ + ll^llxi + 1(^)^)1 sufficiently small, with 

\\T(a,b,e,T,a)\\ Xl <C(\a\ + \\b\\ 2 Xl ), 
\\d a>b T{a,b,e,T)\\ L{Bl)Bl) < C(\a\ + \\b\\ Xl ), 

(4.28) \\d T T(a, b, e, T) \\ L[Bl>Bl) < C(\a\ 2 + \\b\\ 2 Xl ), 
\\d E T{a,b,e,T)\\ L(Bl , Bl) < C(\a\ 2 + ||6||^ + |a|), 

<C(\a\ 2 + \\b\\ 2 Xi + l). 

Proof. Immediate, from Proposition 12.101 and Lemmas I4T31 and 1431 □ 

Proposition 4.7 ( [TZ2] ). Under assumptions (Al)-(A2), (H0)-(H4), (V £ ) 
(alternatively, the weakened assumptions of Remark \3. 5\) . there exists a func- 
tion j3(a,e,T,a), bounded from R A+i to X\ and C 1 from R A+e to B\, with 

(4.29) A 2 (a,p(a,e,T,a),e,T) = 0, 
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(4.30) 



\\{3\\ Xl , ||^|U (K)Bl) <C(|a| 2 + |a|), 
\\d a (3\\L(R^ Bl ) < C\a\, 
||9«/3||lcr 2 ,Bi) < C, 



for | (a, e, T, a)\ sufficiently small. Moreover, for \(a,e,T)\, sufficiently 
small, all solutions of (14.241) lie on the £-parameter manifold {b = /3(a, e, T, a)}. 

Proof. By Corollary 14.41 (14. 24ft is equivalent to the fixed-point problem 



for some a G M. e . By ( I4.28p (i)-(ii). for |(o, e, T, a)\, sufficiently small, T 
preserves a small ball in X\ on which it is contractive in b with respect to 
the weaker norm || • ||^ a . Observing that closed balls in X\ are closed also 
in Bi, we may conclude by the contraction-mapping principle the existence 
of a unique solution (3, which, moreover, inherits the regularity of T in its 
dependence on parameters (a,£,T,a). □ 

Remark 4.8. At the expense of further bookkeeping, we may replace (11.111) 



by \\U\\ Xl := \\e^U\\ H 4, \\d x U\\x 2 := \\dMx, + 1 1 e 2 ^<-> C/ 1 1 ^! and carry one 



further derivative in (a, e, T throughout the analysis, to obtain C 2 depen- 
dence of P(a,e,T,a). Indeed, strengthening (HO) to k = 2r + 1 in (HO), we 
may obtain arbitrary smoothness C r of reduction function (nullcline) /?(•). 

4.3 Proof of the main theorem 

The bifurcation analysis is straightforward now that we have reduced to a 
finite-dimensional problem, the only tricky point being to deal with the £-fold 
multiplicity of solutions (parametrized by a). Define to this end 



with (3 Lipshitz in (a,e,T,a) and C 1 away from a = 0. Solutions (w,v) of 
( 14.51) originating at (a, b) = (a, (3), by (14.121) . remain for < t < T in a cone 



C\ > 0. Likewise, any periodic solution of H4.5[) originating in C, since it 
necessarily satisfyies A 2 = 0, must originate from data (a, b) of this form. 



b = T(a, b, e, T, a) 



(4.31) P(a, e, T, a) := (3 (a, e, T, \a\a), 




(4.33) 



C := {(w,v) : \v\ < C\\w\}, 
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Proof of Theorem ! 1.81 . Defining b = /3(a, e, T, a), and recalling invariance of 
C under flow (14.51) . we may view v(t) multiple 



(4.34) 



v(t) = j(a, e, T, a, t)w(t) 



of w(t), where 7 is bounded, Lipshitz in all arguments, and C 1 away from 
a = 0. Substituting into (I4.5p (i). we obtain a planar ODE 



in approximate Hopf normal form, with nonlinearity N := N w (w,v,e) now 
nonautonomous and depending on the additional parameters (T,a,a), but 
still satisfying the key bounds 



along with planar bifurcation criterion (D £ )(ii). From (14.361) . we find that 
N is C 1 in all arguments, also at a = 0. By standard arguments (see, e.g., 
[HK[ ITZlj ). we thus obtain a classical Hopf bifurcation in the variable w 
with regularity C 1 , yielding existence and uniqueness up to time-translates 
of an ^-parameter family of solutions originating in C, indexed by r and 5 
with r := a± and (without loss of generality) 02 = 0. It remains only to 
establish uniqueness up to spatial translates. 

In the Lax case, we observe, first, that, by dimensional considerations, the 
one-parameter family constructed must agree with the one-parameter family 
of spatial translates. Second, we argue as in [TZ2j that any periodic solution 
has a spatial translate originating in C, yielding uniqueness up to translation 
among all solutions and not only those originating in C; see Proposition 2.17 
and Corollary 2.18 [TZ2J for further details. 

In the overcompressive case, we observe, likewise, that the ^-parameter 
family constructed consists of spatial translates of a smooth (£— l)-parameter 
family of distinct orbits of traveling- wave ODE (II. 4p . then argue as in [TZ2J 
that any periodic solution has a spatial translate lying within a corresponding 
cone C about some member of this [I — l)-parameter family. Constructing 
solutions about each such member by the same technique, we thus obtain 
from their union an (I — 1) parameter family containing translates of all 
periodic solutions on an entire neighborhood of U, and not only in cone C. 
See Proposition 2.22 and Corollary 2.23 [TZ2j for further details. □ 



(4.35) 




(4.36) 



\N\, \d e>TAa N\<C\w 



\d w N\ <C\w 
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Remarks 4.9. The apparent restriction to C 1 regularity caused by factor 
\a\ in the definition of (3 is illusory, since we restrict eventually to the ray 
ci2 = and a\ > 0, on which j3 is as smooth as (3. By Remark I4.8[ the 
latter may be made as smooth as desired by assuming sufficient regularity 
in (HO), and so we can carry out a bifurcation analysis to arbitrary desired 
regularity 

2. Restricting to the "central" solutions a = of the constructed cone 
of solutions yields (3 = 0(\a\ 2 ), and thus the exact Hopf normal form 

(4.37) w = (^fy ^jj w + N w (w, 0, e) + M(w, e, T, t, fit, a), 

where M = 0(\w\ 3 ) and d Wja M = 0(\w\ 2 ). For these central solutions, 
we thus recover all of the standard description of Hopf bifurcation: in par- 
ticular, that e(r) is a pitchfork bifurcation with (de/dr)(0) = (hence 
sgn(d 2 e/dr 2 )(0) is expected to determine stability; see the final paragraph 
of Section [T3]) . 



5 Physical applications 

We conclude by describing applications to various physical systems. To gen- 
erate example systems and waves, it is easier to look for traveling-wave so- 
lutions U(x, t) = U(x — si) with possibly nonzero speed, to be transformed 
to steady solutions by the change of coordinates x — > x — st as described in 
the introduction. Examples satisfying (Al)-(A2), (H0)-(H4) are: 

1. The general Navier-Stokes equations of compressible gas dynamics, 
written in Lagrangian coordinates, appear as 

{v t -u x = 0, 
U t +Px = (iu/v)u x ) x , 
(e + u 2 /2) t + {pu) x = ((k/v)T x + (fi/v)uu x ) x , 

where v > denotes specific volume, u velocity, e > internal energy, 
T = T(v, e) > temperature, p = p{y, e) pressure, and \i > and k > 
are coefficients of viscosity and heat conduction, respectively. For simplicity, 
assume an ideal temperature dependence 

(5.2) T = T(e), 
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independent of v. This assumption can be removed with further effort; see 
Remark I5TT1 Defining U\ := (v), U 2 := (u,e+ \u\ 2 /2), we find that the single 
hyperbolic mode is v , and the associated equation is linear by inspection; 

likewise, the viscosity matrix b = ( V ' V ^ ] is lower triangular with 



k/vc, 

positive real diagonal entries, hence its spectrum is positive real and (by 
Lyapunov's Lemma) there exists A\ 2 symmetric positive definite such that 
A\ 2 b > 0. Conditions (A1)-(A2) and (H1)-(H4) are thus satisfied under the 
mild assumptions of monotone temperature-dependence, 

(5.3) T e > 0, 
and thermodynamic stability of the endstates, 

(5.4) (p v ) ± < 0, (T e )± > 0; 

see, e.g., |MaZ3l IZ2l |Z3] for further discussion. Notably, this allows the 
interesting case of a van der Waals-type equation of state, with p v > for 
some values of v along the connecting profile. 

2. Next, consider the equations of MHD: 
(5.5) 

v t ~ Uxx = 0, 

u lt + {p+ (l/2 y u )(5| + Bl)) x = ((v/v)u lx ) x , 
u 2t - ((l//io)B*B 2 ) x = {{v/v)u 2x ) x , 
u 3t - ((l/fi )B*B 3 ) x = ((u/v)u 3x ) x , 
(vB 2 ) t - {Blu 2 ) x = ((l/afi v)B 2x ) x , 
{vB 3 ) t - (B$u 3 ) x = ((l/an v)B 3x ) x , 
(e + (1/2) (u? + u\ + u\) + {l/2^)v{Bl + Bl)) t 

+ [(p + (l/2/i )(£ 2 2 + Bl)) Ul - (l/^)B{(B 2 u 2 + B 3 u 3 )] x 
= [{p/v)uxU 1x + (fi/v)(u 2 u 2x + U 3 U 3x ) 
+ (k/v)T x + {l/cr^lv){B 2 B 2x + B 3 B 3x )] x , 

where v denotes specific volume, u = {ui,u 2 ,u 3 ) velocity, p = P(v,e) pres- 
sure, B = (B*, B 2 , B 3 ) magnetic induction, B{ constant, e internal energy, 
T = T(v, e) > temperature, and \i > and v > the two coefficients 
of viscosity, k > the coefficient of heat conduction, /x > the magnetic 
permeability, and a > the electrical resistivity. Under assumptions (|5.2|) . 
(15.31) . and (15.41) . conditions (Al)-(A2) are again satisfied, and conditions 
(H1)-(H4) are satisfied (see [MaZ3j ) under the generically satisfied assump- 
tions that the shock be of Lax or overcompressive type, the endstates U± 
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be strictly hyperbolic, and the speed s be nonzero, i.e., the shock move 
with nonzero speed relative to the background fluid velocity, with U\ := (v), 
U2 '■— (u,B,e+\u\ 2 /2 + v\B\ 2 /2fi ). (For gas dynamics, only Lax-type shocks 
and nonzero speeds can occur, and all points U are strictly hyperbolic.) 

3. (MHD with infinite resistivity/permeability) An interesting variation 
of f)5.5p that is of interest in certain astrophysical parameter regimes is the 
limit in which either electrical resistivity a, magnetic permeability fi , or 
both, go to infinity, in which case the righthand sides of the fifth and sixth 
equations of (15. 5p go to zero and there is a three-dimensional set of hyperbolic 
modes (v,vB2,vB 3 ) instead of the usual one. By inspection, the associated 
equations are still linear in the conservative variables. Likewise, (Al)-(A2), 
(H1)-(H4) hold under assumptions (15. 2p . (15.31) . and (15.41) for nonzero speed 
Lax- or overcompressive-type shocks with strictly hyperbolic endstates. 

4. (multi-species gas dynamics or MHD) Another simple example for 
which the hyperbolic modes are vectorial is the case of miscible, multi-species 
flow, neglecting species diffusion, in either gas- or magnetohydrodynamics. In 
this case, the hyperbolic modes consist of k copies of the hyperbolic modes 
for a single species, where k is the number of total species. Again, the 
associated equations are linear, and (Al)-(A2), (H1)-(H4) hold for nonzero 
speed Lax- or overcompressive-type shocks with strictly hyperbolic endstates, 
under assumptions (15. 2p . (15. 3p . and (15. 4p . 

Remark 5.1. Assumption (15.21) may be removed in all cases by working in 
the more general framework described in Remarks 11.31 and 13.51 with variables 
U x = (v), W 2 = (u,T) (resp. W 2 = (u,B,T)). 

A Lagrangian vs. Eulerian formulation 

In this appendix, we discuss a bit further the role of structural assumptions 
(Al)-(A2) in the energy estimates of Section [3J Estimate (13.61) holds by a 
similar, quasilinear version of the argument of Proposition I3~2"} for the general 
class of hyperbolic-parabolic systems considered in [Z3J. However, (13. 7p is 
much more delicate, clS clXG the related estimates (13. 61) . (13. 7p . and (13.141) . 

To see why, note that (13.71) is essentially a variational bound, measuring 
the difference between two solutions, and variational bounds typically cost an 
additional derivative. For example, in the quasilinear symmetric hyperbolic 
case, we have ||^(t)||/f s < C'||^o||h s by standard energy estimates, whereas 
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the same energy estimates applied to the difference between two solutions 
U 1 and U 2 yield, rather, 

WU 1 - U 2 \\ H s{t) < C\\Ul - U 2 \\ H s\\U^\\ HS+ i. 

In the linear case, WU 1 — U 2 \\H s (t) < C\\Uq — Uq \\ 2 Hs follows by superposition, 
and so this issue does not arise; this motivates the assumption of linearity of 
the hyperbolic part of (11.2j) made in (Al). 

Indeed, when (Al) fails, so, typically, does (13.71) . For example, consider 
the Eulerian version (11.91) of the isentropic compressible Navier-Stokes equa- 
tions, rewritten in quasilinear form 



(A.l) 



Pt + Upx + PU X = 0, 

U t + UU X + p~ 1 p'(p)px = VP~ l U a 



Perturbing about the constant solution (p, u) = (1, 0), and assuming without 
loss of generality that p'(l) = 1, we obtain perturbation equations 

/a 2 ) I Pt + up x + u x = -pu X} 

\ U t + p x - VUxx = 0(\p\ 2 + \u\ 2 + \p\\u x \) x 

so lone < C||p||ij-i remains sufficiently small. Applying the same 

energy estimates as in the proof of Proposition 13.21 we easily obtain 

(A.3) \\(p,u)(;T)\\ 2 Ha + [ |K,f)||^i dt < C||(p ^o)||^ 

for s > 1, similarly as in (13. 6p . 
Now, consider the solution p of 

Pt + upx + u x = 0, 

with p, u as determined by flA.2j) . Defining E := p — p, we have 

E t + uE x = pu x . 

Applying the same energy estimates carried out on the variation E in the 
proof of Proposition I3.2[ we thus obtain for s > 1 that 

(A.4) \\p-p\\ H ,(t)<C\\(p ,u )\\ 2 Hs . 
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On the other hand, defining p to be the solution of the linearized equation 

Pt + u x = 0, 

we find 



(A.5) 



\\p(t) - Po\\h° < J \\u\\H°+i(r)dr <t 1/2 (^J \\u\\ 2 H3+ i(r)dr] 

<Ct l / 2 \\(p ,u )\\ H s. 



Likewise, defining p to be the solution of 

Pt + up x = 0, 

and e := p — p we find from similar energy estimates, together with 

^||^e|U 2 = ^||^e||| 2 /2||^e|| L2 , 

that 

(A.6) ||p-p|M*)< / \\u\\H^(r)dr<Ct^ 2 \\{p ,u )\\ H s. 

Jo 

Combining flA.31) . (]A.5p . and (|A.6[) . we obtain 

\\p-p\\ H s(t) = \\p - po\\ H .(t) + O(t 1 / 2 )\\(p ,u )\\ H s, 

which, for t sufficiently small, implies that \\p — p\\H s (t) < C||(Po, M o)||j/s only 
if || p — po||i? s (0 < C 1 1 (po, uo) \\ 2 H s. But, p, as the solution to a simple transport 
equation, is exactly a displacement 

p(x,t) = p (X(x,t)) 

of po, where X satisfies characteristic equation dtX = u(X), and so in the 
absence of derivative bounds on c^po, \\p(t) — Po\\h s is in general no smaller 
than 2||pq||h«. Thus, ||p — p\\H s (t) < C||(p , Wo)||h« in general does not hold 
in the Eulerian formulation. 

The resolution of this apparent paradox is that the Lagrangian formula- 
tion incorporates nonlinear transport effects into the equations via a change 
of independent variables, so that || ■ \\jja compares different solutions along 
their respective particle-paths rather than at a fixed location, thus eliminat- 
ing the principal contribution computed above. 
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Remark A.l. Some such superlinear bound as (13.21) appears necessary 
for bifurcation analysis based on temporal dynamics, at least as usually 
performed based on autonomous linearized equations. This is an obstruc- 
tion to the application of invariant manifold/bifurcation techniques to hy- 
perbolic systems such as the incompressible Euler equations [LiJ. Like- 
wise, for hyperbolic-parabolic systems, it limits us at present to the La- 
grangian formulation and a single spatial dimension. For related discussion, 
see [TIoTIIIio^lHo^ lD]. 
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